ON THE PROPAGATION OF WEAKLY NONLINEAR RANDOM DISPERSIVE 

WAVES 



ANNE-SOPHIE DE SUZZONI AND NIKOLAY TZVETKOV 

Abstract. We study several basic dispersive models with random periodic initial data such that the 
different Fourier modes are independent random variables. Motivated by the vast Physics literature 
on related topics, we then study how much the Fourier modes of the solution at later times remain 
r—t , decorrelated. Our results are sensitive to the resonances associated with the dispersive relation and 

to the particular choice of the initial data. 



1. Introduction 

In this paper we study several basic dispersive models with random periodic initial data such 
that the different Fourier modes are independent random variables. Motivated by the vast Physics 
literature on related topics (see e.g. @), we then study how much the Fourier modes of the solution 
at later times remain decorrelated, and how much the mean values of the amplitudes to the square 
of the Fourier modes vary with time. Our results are sensitive to the resonances associated with 
the dispersive relation and to the particular choice of the initial data. 

All the models we will be interested in can be injected in the following general framework. 
Consider the equation 

(1) (d t + L)u + sJ(u 2 ) — 0, 

posed on the torus T d of dimension d with an initial datum being a random variable that shall be 
described later. In (Q~|), e« 1 since we want to investigate about the effect of a weak non linearity 
over the behaviour of the statistics related to the random initial datum. We suppose that u is real 
valued and L and / are linear maps which are defined as Fourier multipliers by 

Lu{n) - —ia)(n) u(n), Ju(n) - icp(n) u(n), V n e if , 

where"^denotes the Fourier transform on T d and o), <p : Z d \-> R are supposed to be such that 

(2) co(0,n') = <p(0,n') = 0, Vn'eZ", 

with the natural convention in the case d = 1. We suppose that the variable on T d is given by 
x - {x\, . . . , Xd). Then, thanks to the assumption ([2]), we obtain that we can consider solutions of 
(0Q) such that J T u(t, x)dx\ = 0. We also suppose that a>, tp are odd functions. Observe that under 
the last assumption L and J send real valued functions to real valued functions. Set 

D d = {n = (m,..., n d )eZ d \ n x * 0}. 

For s € R, we introduce the Sobolev spaces H s of real functions having zero x\ mean value : 

H s = [u{x) = ^ e mx u n | u n = vlZ, I u{x)dx\ - 0, V |«| 2s |«„| 2 < oo} 

neZ d neV d 

where \n\ = £y \nj\. In this work we shall always make use of these Sobolev spaces H s since they 
are the ones adapted to our models. In all our examples the equation (fT|) is globally well-posed 
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in some H" and thus there will be no difficulty caused by the problem of the existence of the 
dynamics. 

Let us describe the dispersive models which can be written under the form £[]) we will be 
interested in. They all appear in the modelling of long, small amplitude dispersive waves with a 
possible weak transverse perturbation. The first example is the KdV equation 

d,u + d\u + d x {u 2 ) = 

which corresponds to (Q} in the case d - 1 with u)(n) - n 3 and (p(n) = n (with the convention 
x - jci and n = n\ is the case d = 1). The KdV is globally well posed in H s , s > — 1 (see [9], for 
earlier results we refer to Ell5l [T0l '). 

A second example again in the case d = 1 is an alternative of the KdV model, derived by 
Benjamin-Bona-Mahony (BBM equation) which can be written as 

9 9 

d t u + d x u - d t d x u + d x (u ) = 0. 

The BBM equation corresponds to CD) with -u{ri) - ip(n) = n/(l + n 2 ). The BBM equation is 
globally well-posed in H\ s >0 (see |[T1[TT1). 

Our two dimensional models will be the famous Kadomtsev-Petviashvili (KP) equations. In 
fact there are two models according to the impact of the surface tension. The first one is the KP-II 
equation which corresponds to a weak surface tension and can be written as 

d t u + d Xi u + d~ x \d 2 Xl u + d Xl {u 2 ) - 0. 

The KP-II equation corresponds to (HJ in the case d - 2 with io(n\,ri2) = n \ ~ f^/ni if n\ ± 0, 
oj(0, n.2) = and (p{n\,ri2) = n\. The KP-II equation is globally well-posed in H s , s > (see El). 
Finally, the KP-I equation 

d t u + d Xt u — d~ Xx d Xl u + d Xi (u 2 ) = 0. 

corresponds to (OQ) with co(n\,ri2) = n 3 { + n 2 /ni if n\ ± 0, a>(0, n-i) = and ip{yi\,n-i) = n\. The 
KP-I equation is globally well -posed if the data is in H s , s > 2 (see |7] and also (H). 

Next, we describe the random initial data we shall deal with. With T> c [ = {n € D rf \n\ > 0}, let 
(gn) n eD d be a sequence of independent identically distributed complex random variables such that 

E(g n ) = 0, E(\g n \ 2 ) = l 
and such that there exist to positive constants c and C such that for all y e K., 

(3) Eie^^) < C e cyl , £(e rImfe) ) < C e cyl , 

where E is the expectation. We also suppose that the distribution of g n is invariant under the 
multiplication by e' e with 128 + 0[2^]. Note that under these assumptions, E{g 2 ) is equal to 0. 
Further consequencies of this property will be used in the sequel. 

Remark 1.1. A typical example of random variables satisfying our assumptions are the (complex) 
Gaussian random variables, i.e. g„ = -q^(h n + il n ), with h n , l n 6 Af(0, 1). Another example coming 
from the Physics literature is what is known as random phase approximation, that is, g n is written 
gn - XnAn, where Xn is uniformly distributed on S x and A n is a non-negative random variable 
independent from^„, and E(A 2 ) = 1. In all these examples the symmetry assumption on g n holds 
with any angle 6 + 0. In order to ensure (O, we can suppose that the distribution jj. of A n satisfies 

e yr d[i{r) < Ce cyl . 

Jo 

For instance, the last property holds true if jj. is compactly supported. 
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Next, for n e set g_„ = g^. Let A = (A n ) neI) d be a sequence of complex numbers such that 

(4) J] \n\ 2s \A n \ 2 < oo 

neD* 

for some s depending on L and J such that the equation (Q~|) is globally well-posed in H s . Set for 
all n € T>i, A- n = ~A n . Set 

(5) «o(jc) = Ya 8nA " ein ' X 

neV d 

Thanks to our assumption on (A n ) neT) d , we have that uq e H s almost surely. Moreover, it is real 
valued. Let u(s; t, x) be the solution of 

{d t + L)u + sJ(u 2 ) = 
u\t=o = Uq. 

Consider the expansion of u(s; t, x) as a Fourier series, 

u{e; t, x) - u n (s; t)e m x . 

neD d 

Set 5(0 - e~ tL . Then clearly u(0; t, x) = S(t)uo and 

u n (0;t) = e^ n), g n An. 
In particular, thanks to our assumption on the random variables g n , 



(6) E(u m (0;t)u n (0;t)) = 6';;\A n \ 2 , V t e 



Our aim is to expand the quantity E(u m (s; t)u„(s; t)) in s and see how much © survives in the 
nonlinear setting. 

In order to state our result, we introduce the following notations. We set 

A k n l = a)(k) + u(J) - u){ri) 

which corresponds to the pulsation associated to the three waves interaction k + I — > n when 
k + l - n. Next we set : 

Jo 

Here is our main result. 

Theorem 1. Consider (Q]), in the cases of the KP-II, BBM and the KP-I equations, with initial 
data given by ((5]) of typical Sobolev regularity H s , s > 3/8 /or BBM, s > 2 for KP — II and s > 3 
for KP - I. Then 



E(u m (s; t)u n {s; t)) = 6™\A n \ 2 + d" , e 2 G n (A, t) + e 3 R(s; t, m, n) 
where G n (A, t) is given by G„(A, 0) = and 
d t G n (A, t) = 4<p(n) Y Re(-Fi'\-t))((p(n)\A k \ 2 \Ai\ 2 - (p(kMnf\M 2 ~ (KOI^I W) 

k+l=n 

HE(\g n \ 4 ) - 2)(26%>Re{-F q n q {-t))<p 2 {n)\A q \ A - 4Re(-F 2 n n '- n (-t))<p(2nMn)\A n \ 4 ) 

and besides G n (A, t) andR(s; t, m, n) satisfies the following estimates. There exists C > such that 
for every s € (0, 1], every \t\ < —, every m, n, 

\G n (A,t)\ < C\t\ \n\-K'\ \R{s-t,m,n)\ < Cmin(|n|, N)" 1 \t\{\ + \t\) 
in the case of the BBM equation, with B(s) = 2 + 2s if s > 1/2 and B(s) = \ + 4s otherwise, 
\G n (A, 01 < C\t\ \n\- 2s , \R(e; t, m, n)\ < C max(|n|, |m|)|f|(l + |f|) 
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in the case of the KP -II equation, and 

\G n (A, 01 < Ct 2 \n\ 2 ~ 2s , \R(s; t, m, n)\ < C max(|n|, |m|)|?| 3 
in the case of the KP -I equation. 

It is possible that the estimates on the remainder R could be improved. 

It is remarkable that in the case of the BBM equation, if - —r= and E(\g n f) = 2 then 
G n (A,t) = 0. This goes with the fact, due to the first author lfT2l . that the measure (on H l ^ 2 ~) 



induced by 



JILX 



Vl + n 2 

g n Gaussians, is invariant by the flow of BBM. Indeed, this measure is a renormalization of the 



formal measure 



e^du 



and the H l norm of the solution of BBM is conserved by the evolution. In this particular case for 
A n and g n the terms of higher order should also vanish as shows the next proposition. 

Proposition 1.1. Consider the BBM equation. Let A k = y=f an d let g n - -q^(hn + Hn)> with 
h n , l n € Af(0, 1). Then with the notations of Theorem^ 

G n (A, t) = R(s; t, m, n) = 0, V (m, n), V t, We. 

However, in the proof of Theorem [Q the computation of G„ depends only on E(\g n \ ) and 
E(\g n f), which gives a larger framework for almost remaining decorelated initial data. The as- 
sumptions on the random variables that they have large Gaussian deviation estimates is imposed 
in order get the analytic bounds on \R(s; t, m, n)\. 

Remark 1.2. The idea underlying the computation of G n comes from the theory of wave turbu- 
lence and the notion of statistical equilibrium. Indeed, as stochastic laws invariant through the flow 
of one conservative Hamiltonian PDE tend to be quite rare, and to broaden our views on the topic, 
statistical equilibrium is defined as the next best thing, that is, a law whose moments of order 2, 
i.e. the E{\u n \ 2 ) are unchanged by the evolution in time. With F n ((y m ) meT) d, t) = d t G n (( y/y^) m , t), 
if we replace \A m \ 2 by E(\u m \ ) + 0{e 2 ) in the following expression of d t E(\u n \ 2 ), we formally get 
that : 

d t E{\u n \ 2 ) = s 2 d t G n {A,t) + 0(e 3 ) = e 2 F n ((E(\u m \ 2 )) m ,t) + 0(e 3 ) . 
Hence, by neglecting the remainder because of its order in s, we have a closed equation on the 
E(\u m \ 2 ) detecting statistical equilibrium : 

Vt,Vn, d t E(\u n \ 2 ) - F n ((E(\u m \ 2 )) m ,t) = . 

When one takes the weak limit of F n when t goes to oo, only the resonance terms remain. In this 
sense, 

Vn, limF n ((E(\u m \ 2 )) m ,t) = 

t— >oo 

is the kinetic equation corresponding to statistical equilibrium or KZ spectra in the wave turbulence 
theory. Namely, if for instance E(\g n \ 4 ) = 2 then the equation can be written as 

Vn , ^ (ip{n){E{\u k \ 2 )){E{\ Ul \ 2 )) - <p(k)(E(\u n \ 2 ))(E(\u,\ 2 )) - ip{l){E{\u k \ 2 )){E{\u n \ 2 ))) = . 

/t+/=n,A* J =0 

One can see that it depends only on the E(\u n | 2 ), meaning that the solution is invariant through 
dephasing on each wave length. We did not make a serious effort in either finding solutions of this 
equation or of the more general 

(7) Vn W , G n (A, = 
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but we believe that the solutions of the first one, would they exist, would be consistent with the 
KZ spectra and the actual status of the wave turbulence theory. These solutions would act as a 
statistics into which the difference between E(\u n (e; ?)| 2 ) and its initial value would be negligible. 

Remark 1.3. As a matter of a simple observation, inspired by the discussion on the BBM equation, 
we have that if E(\g,,\ 4 ) = 2 then \Ak\ 2 = <p(k)/k\, is a solution of (0. In such a situation the quantity 
E{\u n {s; t)\ 2 ) is the same as its initial value at t = up to a correction of order e 3 , at least for times 
of order 1 . 

Remark 1.4. It also seems that one should be able to get similar results for 



with k > 1, thus approaching the more general law of the solution instead of only the covariances 
between the amplitudes of the Fourier modes. 

Let us observe that the results of Theorem Q] for the KP equation equations also apply for the 
KdV equation, by considering data independent of the transverse variable %2- The result for the 
BBM and KP-II equations is stronger compared to the result for the KP-I equation thanks to the 
absence of resonance interactions. 

The regularity assumptions of Theorem Q] are more restrictive the ones required by the well- 
posedness results quoted above. It is a natural open question whether in Theorem Q] one can 
cover the weaker regularity assumptions of the well-posedness results. We reckon that some new 
phenomenons may occur at low regularities. 

Let us now explain the main ingredients of the proof of Theorem [TJ The first step is to get 
deterministic bounds on the first two Picard iterations. Similarly to the Cauchy problem analysis, 
the presence of resonances plays an important role in the control of the second iteration. We use 
some algebraic cancellations of the average between the first and the second iterations. Similar 
computations appear in the Physics literature. The main novelty in our work is the control on the 
remainder (once one singles out the first two iterations). Here we use an energy method based on 
a conservation law together with the exponential integrability of the first two iterations for times 
of order < e . 

The remaining part of this paper is organized as follows. In the next section we prove Proposi- 
tion ll.ll In the subsequent sections, we prove TheoremQ] 



E(u mi ...u mk u ni ...u nk )) , 



2. 



Proof of Proposition ! 1.11 



Denote by n the measure on H 



induced by the map 




Denote by O(0 the global flow of the BBM on L 2 defined in Q3). Thanks to H2, 



(8) 




Denote by n„ the projection to the rc'th Fourier mode. Then 
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Using © with F(u) = U m (u)U n (u), we get 



E(u m (e; t)u n {s; t)) = H m (u)U n (u)dfi(u) 

Jhv 2 - 

= \ n m (u")n n {u")dp(u>) 

Jo. 

= S^\A n \ 2 . 

This completes the proof of Proposition 11.11 

3. Deterministic estimates for the expansion at order 2 of the solutions 
In this section uq is a deterministic H s function. Consider (Q~|) with data uq. We suppose that 

«0 O) = a n em X > a ~" = ®" ■ 

Let us expand the solution of (0Q) with data uq at order 2 in e. For simplification in the computations, 
let 

v(s;t,x) = S(-t)u(£,t,x) - ^ v„(e, f)e m x 

neV d 

such that v satisfies : 

d t v = -sS{-t)J({S{t)v) 2 ) 

with initial datum uq. Write then 

v(e; t, x) = uo{x) + eb(t, x) + e 2 c(e; t, x) 

and 

v„(e; t) = a n + sb n (t) + £ 2 c n (e; t) 

with 

b(t,x) = - \ S(-T)J((S(T)u ) 2 )dT . 
Jo 

Then 

b n (t) = -i(p(n) ^ a kai F k n ! {t) . 

k+I=n 

The next no resonance lemma plays a key role in out analysis. 



Lemma 3.1. Let k + I = n. Then KP-II and BBM present no resonances, that is for KP-II 
andforBBM 



|A*''| > 3|m*i/i| * 



k l _ nkl(k 2 + P + M + 3) 



*0. 



(1 + n 2 )(l + k 2 ){\ + I 2 ) 

The proof is a straightforward computation. The consequence of this lemma is that the norm of 
b can be bounded independently from t for KP-II and BBM. 

Lemma 3.2. Suppose that for KP-II, the initial datum uq belongs to H s with s > ^, then there 
exists C independent oft and uq such that 

\W)\\ H > < C|M|^ . 

For BBM, suppose that the initial datum is in H s , with s > | and let o~ such that < o~ < 2s - ^, 
then if s < 1, 

\\b(t)\\H<r < Q\U \\ 2 HS , 
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and if s > 5, 

\W)\\h' < Q\uo\\%. . 
Proof. We use the form of b n to give the bound 

\n\ 2s \b n \ 2 < | n | 2i V(«)| 2 2]|a,a„-,a^||^(0l \F?- l (t)\ . 

k,l 

Then, as s > 0, \n\ s < C s (\k\ s + \n — k\ s ), and using the facts that the sum is symmetric in k, n - k 
and l,n — l and that there is no resonances then \F*'"~ k (t)\ < 

\n\ 2s \K\ 2 < C.Mn)\ J] \a k \ \n - k\ s \a n - k \\ai\ \n - l\ s \a n -i\ ^— \ ^ . 
For KP-II, < JU, thus by summing over « and using a Cauchy-Schwartz inequality : 

|A/ I 

\a-k\ \ai\ 

k,i 



I* < 2 ^7 ^7 y In - *|'k^| I" - fl'lo^l 

w ^ 1/11 » 

Hh- < c s \\ UQ \\ H \ y < c,ii«olllr- J E mrr 



2.v 



and as <i = 2, the series converges as long as s > 5. 
In the case of BBM, we have : 



As for s e [-1, 1], we have : 



w _ nkljk 2 + P + M + 3) 
n ~ (1 + n 2 )(l + £ 2 )(1 + I 2 ) 



|yt| i+1 |/| 1_i < k 2 + I 2 < 2(k 2 + l 2 + kl) , 



we conclude that 



, A «l >r N \k\ s+2 \l\ 2 ~ s >nm( „Jk\ S 
lA » l - C l + n 2 l + k 2 l + l 2 - Q(pin)l W 



Hence, 

^ < ci w ■ 

|A*''| 

Let us now bound the norm of ft in terms of the H s norm of uq. Since 

fe„ = -i(f(n) J] a k aiF k n \t) 

k+l=n 

we have that 

2 

|fe„| < |^(n)| 2j l a ^l l^l - 0" 

fe+/=n l^n I 

Using that for cr > 0, |«|°" < C a -(\k\ (T + 1/1°") and by symmetry of the sum over k and / 

InH^I <CV \kf\a k \ \ai\ ^ . 

^ |A*'| 



k+l=n 



We then use the bound on ^j-r to write 



«rift„i < c 2 i^r- j i« fc i iria/i 



k+l=n 
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and therefore 

\\b\\ 2 Ha < c£ £ \ k ^>k\ \n - k\ s \a n . k \\jr s \aj\ \n - j\ s \a n -j\ . 

» k,j 

By reversing the order of the sums and using a Cauchy-Schwartz inequality on the sum over n : 

\\b\\n<r <C\\u \\H°(J]\kr s \a k \) 

tc 

and since 

( V /2 

2>r*M< 2> |2 "~ 4i ' ||Mo||ffi 

k V k , 

and the series converges if s > \ + we get : 

\\b\\H* < Cllwoll^ ■ 

For s > 1 we simply use |A^''| > \(p(n)\ and an argument similar to the one for KP-II yields the 
claimed bound. This completes the proof of Lemma IX2l □ 

Remark 3.1. The arguments we presented here for the KP-II equation relax the assumption s > 1 
to s > 1/2 in 1(131 

For the KP-I equation, there are resonances and hence a much weaker statement holds. 
Lemma 3.3. For KP-I, it appears that for s > 1 

\\b\\ H s-l < C\t\\\uo\\ 2 H s . 
Proof. Use the expression of b to get the bound 

\W)\\h°-i < J \\j((S(T)u ) 2 )dT\\ H1 - t < C s f ||5(r)«ollz.»l|5(r)«ollff < CM \\uq\\ 2 hs 

since as d = 2 and s > 1, the Sobolev embedding H s c L°° holds. □ 

Lemma 3.4. The map c satisfies : 

d tC = -S(-t)j(2S(t)u Q S(t)b + s((S(t)b) 2 + 2S (t)u S (t)c)) + e 2 2S (t)bS (t)c + s 3 (S(t)c) 2 ) . 

Proof. It comes from a combination of the equations satisfied by v and b. □ 

We now would like to prove that c is of order in e but that its order in time depends on the 
cases, whether the equation displays resonances or not. 

Lemma 3.5. For KP equations, one can bound the L 2 norm of c. In the case of KP-I ( with 
resonances), it comes if s > 3 

\\c(t)\\ L 2 < c(t 2 \\u \\ 3 HS +e k| 3 || M oll^)^ £klll " oll " l(1+e| " ll " olbl) - 
And for KP-II, it comes if s > 2 

\\c(t)\\ L 2 < C\t\(\\u \\ 3 HS +e || M oll^> C£klll " oll " l(1+e|l " oll '" ) - 
For BBM, the relevant quantity is the H l norm of c, it comes if s > 3/8 : 

< C\t\(\\u Q \\ 3 HS +£ || Mo ll^)^ e|f|ll " oll " !(1+e|l " oll " s) ■ 
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Proof. Calling E{t) = \\\c{f)\\ L 2 for KP and E(t) = l\\c(t)\\ H i for BBM, 

E(t)d t E(i) = 1 + 11 + 111 

with 

/ = -/«, (2 S < ( «(* + ^«»>), 

and 

For KP equations, it appears that 

7(0 < C\\S(t)c\\ L 2\\d Xl (2S(t)u S(t)b + s(S(t)b) 2 ) || L 2 

and therefore 

7(0 < CE(t)(\\d Xl S(t)u \\ L 2\\S(t)b\\L« + \\d Xl S (t)b\y \\S (t)u \\ L " + \\d Xl S (tMvWS (t)b\\L~) 
Using that the 77 s norms are invariant through the flow 5(0, as s > 1 in both cases, 

\\d Xl S(t)u \\ L 2 < \\S(t)u \\ H i = \\u \\ H i < \\u \\ H * 

7(0 < CE(t)(\\u \\HA\S(t)b\\L« + \\b\\ Hl ||5 (/)moIIl- +s\M# \\S (*)*IIi«0 
and using the fact that J fd Xl (fg) = \ J f 2 d Xl g, for KP we have 

77(0 = - J (S(t)c) 2 d Xl sS(t)u - J (S(t)c) 2 s 2 d Xl S(t)b 
77(0 < Q\S(t)c\\ 2 L2 (s\\d Xl S(t)u \\ L - + £ 2 \\d Xl S(t)b\\L~) 

and thus 

77(0 < CE{t) 2 (s\\u \\ w +s 2 \\d Xl S{t)b\\ L «) . 
Then, for KP-I, use the fact that for s > 3 (s - 2 > 1), 77 v ~ 2 injects itself in L°° 
\\d Xl S(t)b\\L~ < C\\d Xl S(t)b\\H°-i < Q\S(t)b\\H-i ^ C\t\\\u Q \\ 2 HS 

and thus 

8 t E(t) < C(\t\ llHollfp + et 2 \\u \\ 4 HS ) + CE(t)(e\\uo\\ H . + e 2 \t\ \\u \\ 2 HS ) . 
With a Gronwall lemma, 

llc(on L2 < c(^ 2 i| M0 ll^ + £ kl 3 ll"oll 4 ff > £C|?lll " oll " I+£2?2||Ho11 -' . 

For KP-II, use the fact that for s > 2, H s ~ l c L°°, 

\\d Xl S(t)b\\L~ < C\\d Xl S(t)b\\ H *-i < Q\S(t)b\\H> < C\\u \\ 2 H s 

and thus 

8 t E(t) < c(\\u \\ 3 HS + e\\u \\ 4 hs ) + CE(f)(e\\uo\\H> + £ 2 \\u \\ 2 H s) 

and therefore 

\\c(t)\\ L 2 < c(|/|||«oll^ + s\t\ \\u \f H .)e £ ^ H-ollWM ll»oll|. . 
For BBM, c satisfies : 

2(1 - d 2 x )d t c = -S(-t)d x (2S(t)uoS(t)b + s(S(t)b) 2 + e2S(t)u c + s 2 2S (t)bS (t)c + e 3 (S(t)c) 2 ) . 
Since s > |, we can choose a in ]|, 2s - ^[ if s < 1 and cr = s otherwise. We have then that : 
d t \\c(t)\\ Hl < c(||5(0m 5(0^IIl2 + s\\S(t)b\\ 2 L * + \\c{t)\\ H M\S {t)u Q \\ L 2 + £ 2 ||5(0^II L 2)) 
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and as s, cr > 1 /4, the Sobolev embeddings H s c L 4 and H°~ c L 4 hold, 

\\S(t)u S(f)b\\ L 2 < \\S(t)uo\\ L 4\\S(f)b\\ L 4 < C\\S(t)uo\\ H A\S(t)b\\ H , < C\\u \\ 3 HS 

Finally, 

wcmy < c\m\m\\ 3 H . +£ii M oii^)e |f|(e||Holb!+£2||i,oll ^ ) . 

This completes the proof of Lemma [331 □ 

Remark 3.2. One may also establish estimates for higher order derivatives of c by the classical 
energy method. This method does not give the cancellation of the term /// above and thus the 
restriction of the time for which the estimate holds depends on uq and thus on the probability event 
of which uo is a representation. In particular, it is not clear to us how to exploit in general such 
an estimate in the context of the study of the decorrelation of the Fourier modes of the solution. 
Nevertheless, by using random variables g n with values in a compact set, we should be able to 
use the energy method with a time of validity that does not depend on the probability event. For 
instance, one can use g n = XnA n where Xn is uniformly distributed on S 1 and is independent from 
A n , where A n is non-negative, compactly supported and E(A 2 ) = 1. 

4. Probabilistic properties 

In this section uo is given by ([5]). We now want to prove that until time of order e , the maps 
a, b and c are of order in e. For that, we use the following proposition : 

Proposition 4.1. There exist C, c two positive constants such that for allR > 0, the probability for 
the initial datum to have a H s norm bigger than R satisfies : 

P(\\uo\\h> >R)< Ce~ cRl . 
Proof. We first observe that (0 together with the zero mean value assumption imply that 

(9) E{e yKe{8n) ) < e cy \ E(e yJm(gn) ) < e°f . 

First, we notice that thanks to (0, we only need to get (0 for small value of \y\, say \y\ < 1 Next, 
we apply (0 with y - ±a to get 

^orlRefo,)!) + E ( e a\1m(g n )^ < OQ _ 

Now, we use that there exist two positive constants Ci and C2 such that for every \y\ < 1 and every 

x e R, 

\e yx - 1 - yx\ < C l7 2 e CM . 

Thanks to the zero mean value assumption on g n , the above analysis implies that there exists a 
constant A such that 

E( e yRe{g " ] ) < 1 +Ay 2 <e cy \ 

provided c > A. A similar argument applies for the imaginary part of g n . Thus, we indeed have 
© and we are in a position to apply [4, Lemma 3.1]. 

By separating the real and the imaginary parts, using (4l Lemma 3.1], we obtain that there exist 
two positive constants C and c such that for every y > and every sequence (a n ), 

P(\J]a n g n \>y)<Ce- c > 2 ^ n K\ 2 ) 

n 

Indeed, if a n = a„ + if3 n and g n = h n + il„, then, 
I ^ a n gn\ < I Y_i a " hn \ + I ^ 

\YjPnhn\ + \YjPnln\ 
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and therefore, 

P{\Yj a nSn\>y) < P(\Yj a n h n\>y/4) + P(\Yj a n^\>y/^ 

n n n 

+P{\Y u PnK\ > y/4) + PQ £ fi„ln\ > y/4) 

n n 

and then we can apply the [4, Lemma 3.1] on each term of the right handside. Remark that since 
the g„ are independent from each other, so are the h n and the l n , even though h n is not necessarly 
independent from l n . 

We deduce from that that the L q norm (in the probability space ) of ^ & n g n satisfies : 

\\Yj a n8n\\ q Lq <(CqY,^nf) q ' 2 

n n 

with C independent from a n and q. Indeed, this property is due to a change of variable and an 
induction on q. First, we have that : 

\\Y, a nSn\\ q Lq = J qy q - l P{\Y, a nSn\>y)dy< J Cqf^e^'^dy 



With z = 



with 



V2M 2 ' 



\Y, a n8nW%<(Y, \an\ 2 f 2 C(q) 



_2 



C{q) = C J qz q - l e- cr dz 
By integration by parts, we get : 

CO? + 2) = ^C(q) 
2c 

and then using that C(q) is bounded uniformly in q for q e [1, 3], we get 

C{q)<c(£^' 2 <{Cqf 12 

and consequently 

(10) \\Y J a nSn\\ q Lq <(CqY J \ a n\ 1 ) ql2 

n 

Then, we use that : 

P{\\uo\\w >R) = P{\\uo\\ q HS > R q ) < R-«E(\\uo\\%) = ^IImoIIJ,^ 

where L q p denotes the L q norm in the probability space and H s x the H s norm in the physical space. 
For q > 2 and thanks to Minkowski inequality, 

n«oii4,Hi = ii Z i^w^ii^ ^ ii^n^^ii^ 

Hence, using (flOl ) and (0]), we get 

\\m\y p , Hi < n(c?2 w 2 w 2 ' s ') 1/2 i^ = c v^(2 u„i 2 n 2i ) 1/2 < c 

This in turn implies the bound 

P(\\uo\\h* >R)< (|f f /2 ■ 
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Set q(R) - e~ lR ^- such that what inside the parenthesis in the above expression is equal to in 
the particular case q = q(R). If R is such that q(R) > 2 then we have that : 

P(\\m\\H*>R)<e- qiR)l2 = e- cR \ 

Let /?o be defined by 2 = e~ Y jr, i.e. Ro = y2eC. For R € [0,Ro], we can simply write 

P(\\u \\h° > R)e cRl < e cR o 

Therefore 

PQuoWh' >R)< e cR » e~ cRl , VR>0. 
This completes the proof of Proposition 14.11 □ 

As a consequence of Proposition 14.11 we get the following exponential integrability statement. 
Lemma 4.2. There exists So > such that 

We deduce from Lemma 1331 and Lemma [4721 . that as long as \t\ is bounded by e _1 the norm of 
c can be bounded in probability. Indeed, 

Proposition 4.3. For all p > 1 there exists C p such that for all \t\ < we have the bounds : 

E(\\cV)\\ p l2 ) 1/ p < C p t 2 

forKP-I, 

EQmw^yp < c p \t\ 

for KP-II, 

E(]\c(t)\\ P Hl ) 1/p < C p \t\ 

forBBM. 

Proof. This comes from the fact that 

E(\\u \\ p HS ) 

is bounded for all p and that 

£YgC£|f|||«o||HJ\ ^ J7( e c£ 2 \t\\\u \\l s ^ ^ ^gCfi^llMoll^ 

are bounded uniformly in t as long as t ce is less than the 8q defined in Lemma l4~2l □ 
We next collect some properties of the random variables (g n ). 

Lemma 4.4. Under our assumption on (g n ), with the nj belonging to D d , E(g ni g n2 ) = 6„" 2 and 
E(g ni gn 2 gn 3 ) = 0. Moreover E(g ni g m g n3 g n4 ) = 0, unless n\ = -nj for some j e {2,3,4} and 
n k - ~ n l far the two indexes k,l in the set {1,2,3,4}/{1,7'}. Moreover E(g ni g n2 g n3 g n4 ) = 1 if 
n\ t n k andn { ± -n k . Finally E(g ni g n2 g n3 g„ 4 ) = E(\g ni f) if n x = n k orn { = -n k . 

The proof of this lemma follows by using the independence assumption via a careful case by 
case study. In particular, we use that under our assumption of symmetry of the distribution we 
have that 

E(gl) = E(\g n \ 2 g n ) = , E{g A n ) = E{\g n \ 2 g 2 n ) = E(\g n \ 2 g; 2 ) = E(g?) = . 
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5. Expansion of the covariances 

In this section, we complete the proof of TheoremQ] Using the fact that a, b, and c are of order 
in s as long as t < £ _1 we would like to develop the covariances of the amplitudes of the different 
wavelengths. Let d™(t) be defined as 

(%(f) = E(v^(t)v n (t)) . 
Then we have the following statement. 
Proposition 5.1. We have that 

d t d™(t) = 5™£ 2 d t G n (A, t) + £ 3 r( £ ; t, m, n), 
where t r(e; t, m, n) satisfies the bounds for R(s; t, m, n) announced in the statement of Theorem^ 
Proof. Let us compute the time derivative of d%(t) . Since v„ satisfies 

d,v n {t) = -is<p(n) J] v k vie iA »' , 

k+l=n 

it comes 

(11) d t d%(f) = ie<p(m)E( ^ vmvne-^"') - is<p(n)E( £ v k v,v m e iA " J ') . 

k+l=m k+l=n 

In the cases of the KP equations, the term of last order, that is the term of order 7 in s will involve 
three occurrences of c, and since we only have a bound for c in L , we will not be able to bound 

y(n)E[ ^ c k c,c m e iA '' Jt ) 

k+l=n 

by some function depending on the time and not onn,m (see Remark 1X21 ) 

By inserting v n = v„(e) = a n + sb n + s 2 c n {e) in (fTTT ) we distinguish different cases according to 
the power of s. 

First, it is clear that the term of order in the expression of d t d™ is 0. 

Then the term of order 1 is also since it involves three occurrences of a : ataia„, and a^ = A^gk- 
Thus, we can apply Lemma l4~4l This cancellation is frequently used in the Physics literature on 
the subject. 

We will describe the term of order 2 later. 

The term of order 3 involves combinations of 1 c and 2 a or 2 b and 1 a. Hence, in the KP-I 
case it is less than Cmax(|?i|, \m\)t 2 . In the case of KP-II, because of the different estimate on 
b, it is less than Cmax(|?i|, |m|)|f|. A similar analysis applies in the BBM case to get the bound 
C{nao.{\m\,\n\)- l \t\. 

Let us describe this bound in the particular case of combinations between 1 c and 2 a, the 
other ones resulting from similar computations. If we replace one occurence of v by c and two 
occurences of v by a in the expression 

i(p{m)E( Vytv/v ;i e" lA "' f ) - i<p(n)E( J] v t v/V A *') 

k+l=m k+l=n 

we get 



itp(m)E( ^ (a k aic n + a k cia n + c k aia n ) e iA " ; f ) - 

k+I=m 

i<p(n)E( ^ (a k aic m + a k c{d m + c k a{d m ) e' A " ') . 



k+l=n 



14 



ANNE-SOPHIE DE SUZZONI AND NIKOLAY TZVETKOV 



For KP, we can bound the L norm of c, hence, since by a Cauchy-Schwartz inequality on the a 

a k ajc n e' 1 



k+I=m 

and by taking its expectation and circular arguments for the other terms, we get the bound on this 
term: 

3(\<p(n)\ + \<p(m)\)E(\\c\\ 3 L2 ) 1/3 E(\\u \\ 3 Hl ) 2/3 . 

For KP-I, £(||c|| 3 ,,) 1/3 is bounded by Ct 2 and for KP-II it is bounded by C\t\. Hence, as ip(n) = 
this term is bounded by Cmax(|n|, \m\)t 2 for KP-I and Cmax(|n|, \m\)\t\ for KP-II. For BBM, 
£(||c|| 3 2 ) 1/3 < ^(llcll^) 1 / 3 is bounded by C(l + |?|). Hence, as \cp(n)\ < \ri\~ 1 , this term is bounded 

by C(mm(\m\,\n\Y l \t\- 

The third order in e also involves combinations of 2 b and 1 a. In this case, the order in time 
for KP-II and BBM is 0. Hence, the term of third order is bounded by Cmax(|«|, \m\)t 2 for KP-I, 
Cmax(H, |m|)(l + \t\) for KP-II and Cmin(|?i|, + |*|) for BBM. 

The term of order 4 involves combinations of 1 c, 1 b and 1 a or 3 b. Hence, in the KP-I case it 
is less than Cmax(|«|, |m|)|?| 3 . A similar analysis applies in the KP-II and BBM cases. 

The term of order 5 involves combinations of 1 a and 2 c or 2 b and 1 c. Hence, in the KP-I 
case it is less than C max(|n|, |m|)|?| 4 . Again a similar analysis applies in the KP-II and BBM cases. 

The term of order 6 involves combinations of 1 b and 2 c. Hence, in the KP-I case it is less than 
Cmax(|«|, |m|)|/| and a similar analysis applies in the KP-II and BBM cases. 

Finally the term of order 7 involves combinations of 3 c. Hence, it is less than C max(|?i|, |m|)|f| 6 
in the KP-I case and Cmax(|?i|, |m|)|?| 3 in the KP-II case. 

Since t is less than e , we have that all estimates in the KP-I case are 0(max(|«|, \m\)E 3 t 2 ). In 
the KP-II case they are 0(max(|n|, |m|)e 3 (l + |f|)) and in the BBM case 0((rmn(|m|, |n|) _1 e 3 (l + \t\)). 

Let us compute the term of order 2. As it involves 2 a and 1 b, two sums of different nature (and 
their complex conjugate when inverting n and m) appear in it : 

C(*) = i<p(m)E(( ^ auaiKe-^"*) 

k+l=m 

and 

W?(t) - i<p(m)E(( ^ b&iane-^*) 

k+l=m 

which appears twice because of the symmetry between k and /. The term of order 2 is therefore 
equal to 

ceo + vl® + 2(wcco + %m ■ 

By replacing b n (t) by its value 

b„(t) = -hpin) 2_j aja q Fil q {t) 
j+q=n 

we get 

C(0 = <p(n)<p(m) J] 2 E&kWja^e-^Ffit) . 

k+l=m j+q=n 

We now recall that thanks to Lemma 14.41 E(ayiia ja q ) is equal to zero unless we can pair the 
indexes. We can not pair k with I or we will have k = I, m = 0, and m + since m\ ^ by 
hypothesis, but we can pair k with j and / with q or k with q and I with j. In both case we have 
n = m. As long as k + I, we have : 

E(a k aiaja q ) = \A k \ 2 \Ai\ 2 
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otherwise 2k = n and 

E{a k aia jaq ) = 6 2 n k \A k \ 4 E(\g k f) . 
Using our assumptions on the random variables and that e~ lA " r F*''(Y) = -F% l (—t) we get 
O) = 26™<p 2 (n) £ \AM{-F k A-t)) + (EQgnf) ~ 2)8™<p\n)8%{-F™{-t))\A q \ A . 

k+l=n 

Next, 

W?(t) = -<p(m) £ 2 <P(k)E(aja q aia n )e- iA ™ lt Ff(t) • 

k+l=m j+q=k 

Here, we can pair j with / and q with « or j with « and g with / but not j with q. In both case, we 
can do the computation (with changing the indexes): 

n=q=k-j=k+l=m 

As long as / + n, we get: 

E(a k aiajag) = \A n \ 2 \Ai\ 2 

but otherwise 

E(a k aiaja q ) = \A n \ 4 E(\g n \ 4 ) . 
Again, using our assumptions on the random variables, we get 

CXO = ~26%<p(n) ^)K| 2 U,|V /A »V*'-'(f) 



k+I=n 



+(E(\g n \ 4 ) - 2)S:^2n)^{n)e-^ Wf)W 



and since 



we arrive at 



= -25>(«) ^ ^\A\ 2 \M 2 e- iA "\-F*\-t)) 

k+I=n 

+(E(\g n \ 4 ) - 2)6^(2nMn)(-F 2 n "'- n (-t))\A n \ 4 . 

Combining the previous formulae implies the claimed expression for the second order. This com- 
pletes the proof of Proposition [5TT] □ 

Observe that 

dZ(t) = e-^^»E(j£(f)u n (f)) 

and therefore 

E(u^(t)u n (t)) - E(i^(0)K„(0)) = eimn)-oKm)) d m {t) _ d m (Q) 

If m - n, it suffices to employ the fundamental theorem of calculus to the function d" t ' and to apply 
Proposition 15. II If m ± n then one has that d™(0) = and hence one may writre 

e i t (co(n)-^m)) d m {t) _ J m (Q) = _ rf m (())) 

and apply again the fundamental theorem of calculus in combination with Proposition 15.11 

Let us now bound G n (A, t). For KP-I, we use the fact that ip(n) = n\ and \F k ' l (t)\ < \t\. Then, as 
the term involving E(\g n \ 4 ) is included in the sum, 



\d t G n (t)\ < C\t\(\n\ 2 \A k \ 2 \A t \ 2 + \n\\A n \ 2 £ \k\ \A t \ 2 ), 

k+l=n k+l=n 

2 1*1 \Ai\ 2 < (\n\ £ \Ai\ 2 + |/p/| 2 ) < C\n\ \\u \\ 2 HS < C\n\, 



k+I=n 
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\n\ 2s \^\ 2 \U 2 < 2 £ \k\ 2s \M 2 \M 2 < 2max(U / | 2 )|| Mo ||^ < C 

k+l=n k+l=n 

and since |/1„| 2 |ti| 2v is bounded, 

\d t G n (t)\ < C\t\ \n\ 2 ~ 2s , \G n (t)\ < C\t\ 2 \n\ 2 ~ 2s • 
For KP-II, use that ip(n) = n x and \F k J(t)\ < • Hence 

l-M/l 2 ,,,2V-I^I 2 



|a , G „ (()|sc(l „ llE ^L + w2E ^), 



,|2 _ |i,|2 



v i^i 2 . 

and since |i„| 2 |«| 2v is bounded : 

\d t G n (t)\ < C\n\- 2s , \G n (t)\ < C\t\ \n\~ 2s . 
For BBM, we have that <p{ri) = so \ip(n)\ < \ti\~ 1 and 

\F n ' (01 < — -r < ^H—ry < Qn\ — 

for all a between - 1 and 1 . 

Let us now compute the bound for G n (t). Set fi(s) = 1 + As if s e]|, 5] and =2 + 25' 
otherwise. We have 

<9,G n (o| < C|<p(n)| J] |Re(-F„ w (0)|(|^(«)||^| 2 U,| 2 + |<p(*)||/U 2 M 2 + |^(/)P*| 2 U, t | 2 ) 

k+l=n 

\d t G n (t)\ < C(I + II) 

with 

/ = |<p(n)| £ |Re(-F«(/)V(n)| |^| 2 U/| 2 

and 

// = \<p(n)\ J] |Re(-F£'(r)V(*)| U„| 2 U,| 2 . 
By symmetry over k and Z, we get that : 

,2 



|„P|/| < C J] l^(n)ll^ ( ' 5K1 |Re(-^'(0)||^| 2 U/| 2 

k+l=n 

Hence, by using the bound on F*' 1 

\nf^\I\ <C J] I^S^lW 

k+l=n W 

For s < \, choose a = -2s, ie fi(s) - 1 + a = 2s, 

\nf s \l\ < C £ |£| 2 ' v |/| 2 'Ul 2 U/| 2 



k+l=n 

-\2s\ j |2 



which is bounded by the sum of the \k\ \Ak\ to the square. For s > 1/2, choose a = -1, 

\nf s) \I\ <CJ| \k\ 2s \l\ 2 \A k \ 2 \M 2 - 



k+l=n 
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For the second term, we have : 

\nf (s) \II\ < Il.i + II.ii 

with, as |«| 2,s |/l„| 2 is bounded 



k+l=n 

and 



Il.i = C\<p{n)\ 2 \kf is) - 2s \Rz{-F k M)\m)U 



2 

I\ 



II. U = C\<p{ri)\ J] \lf^- 2s \Re{-F k J{m\<f(k)U 

We get 



k+l=n 



k+I=n 



With or - 2s for s < 1/2, 



and with a = 1 otherwise, 



We also have 



Il.i <CJ| l/p'V,! 2 < oo 



k+l=n 



Il.i < C £ |/| 2 |i,| 2 < oo 



k+I=n 



ii.u<c ^ ifl* ,)-2 '~ or |jfcr 1 |/y 2 



k+I=n 



with a = 1, we get in the case s < 1/2 that /i(s) - 2s — a = 1 + 4s - 2s - 1 = 2s and in the case 
s > 1/2 that /?0) - 2s - a = 2 + 2s - 2s - 1 = 1 < 2s, hence 



2 

< oo . 



//.»<C£|Z| 2 U| 
Therefore, in the case of BBM, 

< Qt\ . 

This completes the proof of Theorem Q] 

Remark 5.1. In the case of the KP equations, formally, if A% do not depend on k and g n are 
standard complex gaussians then d t G„(A, t) = 0. This goes with the fact that the measure induced 
by 



should be formally invariant through the flow of KP, as it is a renormalization of the formal object 



and the L 2 norm of the solution of KP does not depend on time. However the support of these 
measures in the case of the KP equation contains functions which are too singular for the available 
well-posedness theory. In the case of the KdV equation the measures is supported by H' l ^ 2 ~ . This 
could be a motivation to try to lower the regularity assumption in our approach in the KdV case. 
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